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A new measurement of the bb production cross section in 920 GeV proton-nucleus collisions is 
presented by the HERA-B collaboration. The bb production is tagged via inclusive bottom quark 
decays into J/ip mesons, by exploiting the longitudinal separation of J/ip —> Ul~ decay vertices from 
the primary proton-nucleus interaction point. Both e"'"e“ and j-L~ channels are reconstructed, for 
a total of 83 ± 12 inclusive b —> J/ipX events found. The combined analysis yields a bb to prompt 

J/ip cross section ratio of ~ 0.032 ± 0.005stat ± O.OOlsys measured in the xf acceptance 

(—0.35 < Xf < 0.15), extrapolated to a{bb) — 14.9 ± 2.2atat ± 2.4sys nb/nucleon in the total phase 
space. 

PACS numbers: 

13.85. Ni Inclusive production with identified hadrons 

13.85. Qk Inclusive production with identified leptons, photons, or other non-hadronic particles 
13.20.He Decays of bottom mesons 

24.85. +P Quarks, gluons, and QCD in nuclei and nuclear processes 
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I. INTRODUCTION 

Recent improvements in the theoretical description of heavy flavor (cc, bb) hadroproduction, inspired by 
the availability of increasing amount of data in various kinematic regimes, have provided new insight into the 
physics of hadron-hadron collisions dll Investigations in this field are further stimulated by the current 
and the next generation of heavy ion colliders (RHIC, LHC) where signatures of Quark Gluon Plasma (QGP) 
can be identified only after hadron-hadron collisions in a non-QGP regime are understood with sufficient 
accuracy. 

In the held of fixed target bottom production, several experiments have been performed over the past 
years at Fermilab, GERN and DESY with both pion and proton beams (see H, for recent reviews). 
Several theoretical calculations have become available [ 3 , although the accuracy of the predictions is 
usually quite poor. The large experimental and theoretical uncertainties provide a strong motivation for a 
more precise and accurate measurement of near-threshold bb production, where the sensitivity to theoretical 
models is large and the kinematic region is complementary to that of cc or bb production at very high 
energies at collider experiments. In fact, in the near-threshold region, perturbative QCD can be applied 
and full next-to-leading-order calculations are available, although higher order terms are important and not 
yet fully evaluated. In this regime, measurements have been published by only three experiments: E789 
E771 0 and HERA-B 0. ’ _ 

The HERA-B experiment can measure the bb cross section via the inclusive b —> J/ipX decay mode. 
Interactions are produced by 920 GeV protons in the halo of the HERA beam impinging on target wires of 
different materials at ^/s = 41.6 GeV proton-nucleon center-of-mass energy. The bb production cross section 
(tjA) on a nucleus of atomic number A is obtained via the inclusive reaction 

pA^bbX with bb ^ J/'ipY ^ {e^ e~ / pL~)Y. 

The b hadrons decaying into J/4’ {“b ^ J/4'^ in the following) are distinguished from the large background 
of J/V’ mesons produced directly on the target (“prompt J/'0” in the following) by exploiting the b lifetime 
in a detached vertex analysis. HERA-B is capable of detecting both dilepton decay channels of the J/4^ 
providing an increase of the statistical significance of the measurement and internal consistency checks. 

A first measurement based on limited statistics was published in |ir| . This work presents an improved 
measurement based on larger statistics collected during the 2002-2003 data taking. The larger sample of 
bb events allows us to perform further cross checks, such as a lifetime measurement, a search for specific 
decay channels of B mesons with a J/4 in the final state and a search for inclusive decays of the form 
b J/ij} X. 

The present paper is structured as follows. In Sect. im the basic concepts of the a{bb) measurement 
are discussed. An overview of the detector, the trigger performance and the data samples is given in 
Sect, mil followed by a description of the Monte Carlo simulation (Sect. 113. The methods used for the 
J/4 reconstruction and selection are presented in Sect. M and discussed separately for the muon (J/'f/) —> 
pAp~) and the electron {J/4 e^e~) channel. Sect. IVT| describes the procedure for the identification of 

secondary vertices from b —*■ J /4 decays and their separation from the prompt J /4 background. In Sect. IVlH 
independent confirmations of the b content of the selected events are provided. In Sect. IVIHl the muon and 
electron data are combined into the final result, which is then compared with the existing experimental 
data on a(bb) and recent theoretical predictions. The main systematic uncertainties of the measurement are 
discussed. Sect. HXI provides a brief summary of the analysis. 


II. MEASUREMENT METHOD 

The total bb production cross section in proton-nucleon (pN) interactions, integrated over the complete 
phase space, is the quantity of interest for characterizing the b production rate since it is most readily 
compared to theoretical calculations and other experimental results. However, in fixed-target experiments, 
the proton-nucleus (pA) cross section is measured for the actual detector acceptance (^crA) and then 1) scaled 
to the proton-nucleon cross section taking into account the known (or assumed) nuclear dependence and 2) 
extrapolated to the full phase space by using theoretical models or information from other experiments. 
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In the present analysis, systematic uncertainties related to detector and trigger efficiencies are minimized 
by measuring the bb cross section relative to the J/ij} cross section in the kinematic region in our acceptance. 
The analysis consists essentially of two steps: 1) the J/'0 selection to determine the number (np) of prompt 
J/ij) mesons (Sect. 0 which are copiously produced on the target; 2) the selection and counting of J/V’ 
mesons detached from the primary interaction vertex (Sect. ED, originating from b ^ J/ip decays (Uf^j). 
Following this strategy, the ratio ofbb to prompt J/ip production cross sections in our acceptance {Aa^/Aap) 
can be written as 


AcTp 


“■bb 

np 


ep-e^- -Bvibb^ J/^Xy 


( 1 ) 


with ep = where and Sp'^ are the detection efficiencies (including trigger, reconstruction 

and selection) for J/tj) from b decays and prompt J/ip, respectively, is the detached vertex selection 
efficiency. Br(&6—*■ J jyx) = 2 • (1.16 ± 0.10)% is the inclusive bb —*■ J/yX branching ratio, assumed to 
be the same in hadroproduction as that measured in Z decays. 

All the quantities entering into Eq. JD (except the branching ratio) are evaluated for J/y mesons produced 
in the kinematic range covered by HERA-B. In terms of Feynman-a; {xp) and transverse momentum (pt), 
the kinematic range covers the intervals —0.35 < xp < 0.15 and 0 < pt < 6 GeV/c. This will be indicated 
in the following as the HERA-B acceptance range. 

The prompt J/y cross section (cTp) has been measured at various energies (yi G [10, 200] GeV) and using 
various target materials {A G [1, 197]) both in fixed target and collider experiments. Cp is usually expressed 
by parameterizing the atomic weight (A) dependence of the cross section as Cp = ■ A“, where 

is the J/y production cross section for proton-nucleon interactions and a has been measured with high 
statistical precision as a function of a:;’ and pp by the E866 collaboration m providing an average value of 
a = 0.96 ± 0.01 in our xp range. 

The actual value of crj/^ is needed to obtain the absolute a{bb) cross section and it is known P| that the 
values at energies around ^/s = 20 — 50 GeV scatter strongly. To determine the value of ujj^ at the HERA-B 
energy, a global analysis has been performed on all published J/y cross section measurements^. The 
best value, obtained from a fit on <Jj/.,p{y/s) data with the help of a NRQCD model evaluated at the next-to- 
leading order is crj/^(41.6 GeV) = (502 ±44) nb/nucleon. The main uncertainty is systematic and has been 
evaluated by changing the selection of data submitted to the fit and by changing the free parameters in the 
model. Taking into account that the fraction of prompt J/y mesons produced in the HERA-B acceptance 
range is fp = (83±1)% (see Sect. nvii . we will assume a reference prompt J/y cross section in our acceptance 
which is given by: Actj/^ = fp crj/tp = (417 ± 37) nb/nucleon. 

The A dependence of the crA cross section has not been measured. In this paper, we assume the form 
(jA = a{bb) ■ A"*", with Of, = 1 since no nuclear suppression is expected in open bottom or open charm 
hadroproduction M This expectation is confirmed in H-meson production data M- The same assumption 
is made in all published fixed-target a(bb) measurements. 

The different A dependences of prompt J/y and bb cross sections complicate the evaluation of the cross 
section ratio when combining data sets from different target materials and several approaches are possible. 
Since in the present case suffers from low statistics in some sub-samples, we express the bb to J/y cross 
section ratio per nucleon in our acceptance range as: 


Ract 


A(t(56) 


-bb 


Br(66^ J/yX) ■ S,n*p • ep,i ■ eA- • ’ 


bb,; 


( 2 ) 


where is the total number of detached J/y events found in the acceptance, rip, ep,i and are 
respectively the numbers of prompt J/'0, the J/ip selection efficiency ratios and the detached vertex selection 


^ The global fit includes our own measurement of the J/tp cross section (663di74di46 nb/nucleon Il4l i. We prefer to normalize 
the present measurement to the value given by the global fit since the fit provides a complete summary of all available 
measurements. 
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FIG. 1: Plan view of the HERA-B detector. 


efficiencies measured on the sub-sample i taken with a target of atomic weight Ai. The quantity defined by 
Eq. (|2J) has the advantage of minimizing the dependence on theoretical models and on extrapolations. We 
also quote the bb cross section for full phase space, evaluated as a{bb) = R^a ■ where is the 

fraction of J/V’ mesons from b decays in the HERA-B acceptance range. The quantity is determined 
through theoretical models as discussed in Sect. ITVI 


III. DETECTOR, TRIGGER AND DATA SAMPLE 

HERA-B [I3I3 is a fixed target experiment at the HERA storage ring at DESY. The spectrometer and 
the trigger system were designed for efficient real-time filtering of b hadrons based on an online reconstruction 
oi Jftp mesons. 


A. Detector 

The spectrometer (Fig.^ has a forward geometry, covering from 15 to 220 mrad in the bending plane, and 
from 15 to 160 mrad vertically. The detector features high resolution tracking and vertex reconstruction, 
good particle identification over a large momentum range, and a multilevel dilepton trigger. 

The wire target M consists of two independent stations containing 4 wires each and separated by 40 mm 
along the beam direction. The wires are made of various materials (C, Al, Ti, Pd, W) with dimensions of 
50-100 /im perpendicular to the beam and 50-500 /im along the beam. Up to eight target wires can be 
operated simultaneously and each wire can be moved independently into the proton beam halo to adjust the 
interaction rate. During the data taking, the target operated at rates of about 7 MHz, both in single and 
double wire configurations. 

The Vertex Detector System (VDS) located downstream of the target, consists of 8 planar stations 
with a total of 64 silicon microstrip detectors (~ 50/im readout pitch). It provides high spatial resolution for 
the reconstruction of primary and secondary vertices. A resolution of about 450 /im on Jftp decay vertices 
along the direction of flight was achieved, fulfilling design specifications. 

The main tracking system is located behind a 2.13 T-m bending magnet and extends to 13 m downstream 
of the interaction region. The granularity of the tracking system was adapted to the increasing particle 
density with decreasing distance to the beam. The Inner TRacker (ITR) covers the range u p to 20 mrad 
with Micro Strip Gas Chambers (MSGC) equipped with GEM foils. The Outer TRacker (OTR) |22l| . made of 
honeycomb drift chambers, reached a hit resolution of about 350 fim. The overall tracking system, consisting 
of VDS, OTR and ITR, reached a muon momentum resolution of Ap/p{%) = (1.6I±0.02)0(0.005I±0.0006)p 
where the momentum p is in GeV/c and in the range 5 < p < 80 GeV/c. 
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Particle identification is performed by a Ring Imaging Cherenkov detector (RICH), an Electromagnetic 
CALorimeter (ECAL) and a MUON system. 

The RICH detector provides pion-kaon-proton separation. It uses a 2.5 m long C 4 F 10 radiator volume 
for Cherenkov light emission by charged particles, which after reflection in spherical and planar mirrors, is 
directed to one of two planes of photomultiplier tubes. 

The ECAL is optimized for good electron/gamma energy resolution and electron-hadron discrimi¬ 
nation. It is composed of 5956 independent calorimeter modules built with shashlik technology. The cell 
size increases with the distance from the beam. The detector is instrumented with a fast digital read¬ 
out and a pretrigger system which provides reconstructed clusters with large transverse energy (pretrig¬ 
ger seeds). The energy resolutions for the three ECAL regions (Inner, Middle, Outer) are, respectively, 
oeIE = 20.5%/VE © 1.2%, ueIE = 11.8%/v^ © 1.4% and = 10.8%//E © 1.0% |2|. Spatial 

resolutions, determined with decays, are ranging from 1 to 8 mm depending on the ECAL region and the 
particle energy. 

The MUON system consists of four tracking stations located in the most downstream part of the 
detector at various depths in an iron and concrete muon-absorber. The active elements are mainly con¬ 
ventional tube chambers arranged in double layers. Pad cathodes coupled to the tubes of the two most 
downstream tracking stations are used for pre-triggering purposes. 


B. Trigger Configuration 

The dilepton trigger is initiated by pretrigger (PT) signals from either the MUON system or the 
ECAL [^. Dual pad coincidences in the MUON system or complete ECAL clusters with a transverse 
energy above 1 GeV constitute the basic trigger seeds. These are used by the First Level Trigger (FLT) 
to define a sequence of regions of interest inside the OTR, which are used to initiate a search for track 
candidates originating from^U'0 decays. Track parameters as well as the PT seeds are then sent to the 
Second Level Trigger (SLT) |^, which is a highly configurable software filter that extrapolates found tracks 
through the magnet and VDS, finally applying a mild vertex constraint. 

More specifically, for the data presented here, the FLT required the presence of at least two PT seeds 
and at least one reconstructed track originating from one of the PT seeds. The SLT searched for tracks 
starting from the PT seeds, requiring that at least two complete tracks, with segments in the OTR and the 
VDS, be found and that the two tracks be consistent with a common vertex hypothesis. Moreover, the SLT 
track-seeding algorithm imposed a target constraint by including in the initial track parameters a point near 
the center of the active target wires with uncertainties large enough to accommodate tracks from all targets 
and also to allow for the finite B lifetime. Accepted events were read out and sent to an online reconstruction 
farm [^. 

The J14’ trigger worked constantly at rejection factors around 4 x 10^, matching the typical data recording 
rate of 120 Hz (to tape). The experiment routinely achieved an event yield of about 1200 J/4 mesons per 
hour at an interaction rate of about 7 MHz. 


C. Data Sample 

The data sample used in this analysis was collected between October 2002 and March 2003. During this 
period, 164 million dilepton-triggered events were recorded containing about 300,000 J/ip mesons. These 
are distributed almost equally between the dimuon and dielectron channels. 

Data were taken in nine different wire configurations of single or double wire runs, for a total of 14 
equivalent single-wire samples. The wire materials used were carbon (A= 12 , «64% of the J/4 sample), 
tungsten (A=184, «27%), and titanium (A=48, «9%). 

During the data taking period, the detector was operated under constant monitoring. The data quality 
was assessed both online and offline. Only runs with good performance of all the main detector components 
and stable trigger conditions are used. Five periods with constant experimental conditions were defined and 
the data were grouped accordingly. 
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IV. MONTE CARLO SIMULATION 


Monte Carlo (MC) simulations are used to determine the efficiency terms entering into the cross section 
equation (Eq. (21), to determine the criteria to select the b —>■ J/ip candidates, and to understand the nature 
of the background. The simulation includes our best knowledge of the physics of the processes under 
investigation, of the detector, of the trigger and of the event reconstruction. 

To best describe interactions in the HERA-B environment, the Monte Carlo generator exploits a combi¬ 
nation of features of different standard tools. The simulation of heavy quark (Q) production is obtained by 
generating the basic process pN QQX with Pythia 5.7 [s^ and hadronizing the heavy quarks using the 
Jetset 7.4 1^ package. In the second step, the remaining energy and momentum of the collision is used 
by Fritiof 7.02 to provide the underlying inelastic event and to generate further interactions inside 
the target nucleus. After these steps, the generated particles are passed through the Geant 3.21 package 
[s^ for detector simulation. 

For an accurate description of the kinematic characteristics of prompt J/'i/' production, experimental data 


on prompt J/^/> differential cross sections da/d(p\) and da/dxp are used to tune the standard generator: 
all the produced events are weighted to match the available proton-silicon measurements on differential J/tp 
cross sections. Since these are available in the positive xp region only, we assume a symmetric xp distribution 
of prompt J/ijj production: da/dxp = /3-(l— where c = 6.38±0.24 |^. This extension conforms both 
to all the basic charmonium production models and to our own data on the J/ip differential cross sections 
|3(i|| . The model dependence of the generated pp spectrum is of minor importance since the acceptance 
for the J/tp pt is essentially flat. The influence of a possible J/ip polarization within the limits given by 
experimental results |35L l37l| is taken into account in the systematic uncertainties of the prompt J/p} MC. 
With this model, the fraction fp of prompt J/ip mesons produced in the HERA-B acceptance range is 
A = (83 ± 1)%, where the error is due to the uncertainty in the exponent c of the proton-silicon data fit 


Since, for open b production, there are no measured differential cross sections, we use a bb production 
model which is mostly based on NRQCD theory . As in the case of charmonium production, the Pythia 
generated bb events are weighted with a model containing several contributions: 1) the generated b quark 
kinematics {xp and pp), as given by the computation of Mangano et al., using the most recent next-to- 
leading-order MRST parton distribution function with a b quark mass of mb = 4.75 GeV/c^ and a QCD 
renormalization scale p = + P^] 2) the intrinsic transverse momenta of the colliding quarks, smeared 

with a Gaussian distribution resulting in (fc^) = 0.5 GeV^/cM^; and 3) the b fragmentation function, 
described by a Peterson shape a parameter e = 0.006 j42l|. The subsequent b hadron formation and 

decay are described by the Pythia default parameters. Possible b hadron polarization and J/ip polarization 
in b decays are assumed to be null. With these assumptions, the fraction of detached J/ip mesons produced 
in the HERA-B acceptance range (see Sect. nH) is fbb = (90.6 ± 0.5)%, where the error is determined in the 
systematic studies described below. 

The sensitivity of the measured cross section ratio Ra<t to the bb MC model parameters has been deter¬ 
mined by varying: the parton distribution functions (from MRST to CTEQ5 [43), the b quark mass (in the 
range mt, G [4.5,5.0] GeV/c^), the QCD renormalization scale (from 0.5^/ml + p^ to 2^/ml -bpy), the frag¬ 
mentation function (from the Peterson form with parameter e G [0.002, 0.008], to the Kartvelishvili 

form |45| with parameter = 13.7 ± 1.3 |44|1. the intrinsic transverse momentum distribution (with (fc|.) 
in the range [0.125, 2.0] GeV^/c^) and the fraction of b baryons produced in the b hadronization process in 
the range [0,12]%. The observed variations in the detection efficiencies have been included in the systematic 
uncertainty. No significant dependence (< 1.5%) of the cross section ratio has been found on the momentum 
of the J/ip mesons in the b hadron rest frame nor on the polarization of J/ip mesons coming from b decays. 

The Geant tracking step is performed with the help of a detailed description of the detector which includes 
both the active (instrumented) and inactive (support structure) elements. The detector response is simulated 
by reproducing the digitization of electronic signals and of the read-out chain. A status (masked/dead/alive) 
and an efficiency is associated to each channel to take into account the measured efficiency in data as well 
as the faulty/noisy/masked channels. This information was produced for each data-taking sub-period (see 
sect. rrrm . The MC events are subjected to a full trigger simulation and reconstructed with the same code 
as the data. For each data sub-sample, MC samples were produced for the signal channels and background 
sources to evaluate the efficiency terms in Eq. © for the specific running conditions. 





V. J/V> EVENT SELECTION 


The first stage of the analysis is devoted to the determination of the number of J/ip mesons produced 
at the target, which is dominated by the prompt production. This is done in several steps: selection of 
well-measured tracks, lepton identihcation, lepton pair selection and J/ip candidate counting. Most of the 
selection criteria are common for the muon and the electron final states and are described first. Other details 
specific to the different lepton types (for example, particle identihcation and counting of prompt candidates) 
are discussed separately. 

The cuts applied for prompt J/ip signal selection mostly rehect the detector acceptance and requirements 
already explicitly or implicitly imposed by the trigger algorithms. Therefore the lepton track search initially 
uses information provided by the trigger. The parameters of the SLT track candidates that resulted in a 
trigger decision are given as seeds to the offline reconstruction program. If a seed successfully results in 
a reconstructed track, it is flagged as a ‘trigger track’. Only these are considered in the search for J/ij} 
candidates. 

Further selection is performed on the offline reconstructed tracks, which proht from improved reconstruc¬ 
tion due to more accurate calibration and alignment constants compared to those used online. Basic cuts 
on the number of hits in the VDS and in the main tracker (OTR-I-ITR), and on the probability of the 
track fit are used to select only well-reconstructed tracks and to reject mismeasured and ghost tracks. Weak 
cuts on momentum (6.0 < p < 200 GeV/c) and transverse momentum (0.7 < pr < 5.0 GeV/c) are applied. 
The lower cuts mostly reflect the detector acceptance and the trigger requirements. In the muon case, they 
also improve the rejection of hadrons traversing the absorber and of muons from tt jK decays in flight. The 
upper cuts discard tracks with momenta above the expected range for Jjtj) decay products. 

To complete the single lepton selection, requirements on the relevant particle identihcation specihc to 
muon or electron channels are imposed and are discussed separately in Sect. EH and EE 

When two lepton candidates with opposite charge are selected, a vertex ht is performed. Only dilepton 
candidates with a probability greater than 1% are considered further. This requirement ensures that the 
tracks originate from the same vertex and helps to reject combinatorial background and background from cc 
and bb double semileptonic decays in the subsequent detached vertex analysis. Moreover, the dilepton is 
required to be between —0.35 and 0.15 to match the detector acceptance. Only one dilepton candidate per 
event is accepted. When more than one dilepton candidate survives all cuts, the one with the best particle 
identihcation for both leptons is chosen. 


A. J/tp p 

The identihcation of muon tracks is simplihed by the fact that only muons have a signihcant probability 
of penetrating through the absorbers of the MUON detector. Therefore a minimal set of cuts on signals 
from this detector (hit multiplicity and muon likelihood C^) is sufficient for robust muon identihcation. In 
case more than one muon pair is selected, the pair with the highest product of muon likelihoods is chosen. 

The dimuon mass spectrum obtained with this selection is shown in Fig. [2^. The number of prompt J/tp 
candidates is determined by htting the mass distribution to a function derived mainly from Monte Garlo 
studies ^ 5 . This function includes: 1) a description of the non-Gaussian shape of the J/ip signaP, 2) a 
radiative tail due to the J/ip —> pL~^ decay, 3) a Gaussian shape for the 'ip{2S) signal and 4) an exponential 
function for the background. 

The mean value of the J/ip mass is 3093.5 ±0.2stat MeV/c^, and the full width at half maximum (FWHM) 
of the background subtracted J/'ip peak is 90 MeV/c^. The shift observed in the J/ip mass with respect 
to the PDG value 0 is due to small systematic effects in the tracking and in the internal alignment of 
the detector. The total number of reconstructed J/ip mesons is np{p~^p~) = 148200 ± 500. This number 
depends slightly on the fitting function and mass interval considered. These uncertainties are included in 
the systematic error ISect. IVimi . 


2 


A symmetric function composed of three Gaussians is used, to take properly into account the details of the momentum 
resolution in a complex apparatus such as HERA-B. 
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B. J/i/> ^ e+e- 

The electron and positron selection is affected by larger background contributions than the muon case, 
mostly due to pions interacting in the ECAL and hadrons overlapping with energetic neutral showers. 
Requirements on the lepton identification are therefore more demanding than in the muon case. Additional 
corrections are used to improve the accuracy of the electron momentum estimate. 

The electron track candidates are selected by requiring an ECAL cluster with a transverse energy greater 
than 1.0 GeV, matched to the track, both at the trigger and the offline analysis levels. In the non-bending 
direction (j/), the distance Ay between the track position extrapolated to the ECAL and the reconstructed 
ECAL cluster must satisfy the requirement \ Ay/aAy\ < 3.0, where a Ay is the resolution on Ay for clusters and 
tracks from e^. The typical resolutions are in the range 3—10 mm depending on the ECAL region. Cluster- 
track matches from random combinations and from hadron interactions in the calorimeter are substantially 
suppressed by this cut since their Ay distributions are wider by more than a factor two compared to 
tracks. 

For a precise measurement of the electron momentum, some corrections to the value provided by the 
tracking system are needed. During their passage through the material in front of the ECAL, electrons can 
emit bremsstrahlung (BR) photons and lose energy due to ionization processes. ECAL clusters due to BR 
photons emitted before the magnet are searched for in the ECAL at a location given by the extrapolation of 
the initial track direction to the ECAL position. If a BR cluster with a measured energy Ebr > 1 GeV is 
found, the energy of the cluster is added to the electron momentum at the production point. The correction 
(« 1 — 2%) for further energy losses is evaluated stochastically via Monte Carlo simulations. 

Electron identification in the ECAL is mainly based on the ratio of the cluster energy (A) to the momentum 
obtained from the tracking (p). Typical E/p distributions for signals have a Gaussian shape with mean 
0.99 and width a ~ 6%. Since the requirement on E/p is also crucial in the selection of the detached 
signal, it was included in the optimization procedure described in Sect. El The obtained cut values are 
—2.0 a < Ejp— 0.99 < 3.5 a. If more than one electron-positron pair is found in one event, the combination 
with the best E/p ratio is taken. 

Requiring that at least one of the electrons from a J/ip candidate has an associated bremsstrahlung 
photon (BR tag) or that both leptons have BR tags increases the signal-to-background ratio from 0.6 to 
3 and 9, respectively, providing a powerful method to isolate clean samples. These are used for studies of 
resolutions (for example, of Ay and E/p) and for cross checks. As explained in ^3) the probability that a 
BR photon is emitted and detected in a separated ECAL cluster (bremsstrahlung tag probability) can be 
measured redundantly by counting J/i/) mesons for different tag criteria. For the current setup, the average 
tag efficiency (ebr) is Ebr = 0.30±0.01stat±0.0Isys per track. The same analysis, performed on MC events, 
yields a MC tag probability of e^^ = 0.29 ± O.OIstat ± O.Olgys, in agreement with data. This gives further 
confidence in the quality of our detector description and MC simulation, in particular in the crucial region 
in front of the magnet, where the triggered electrons cross about 0.07 radiation lengths of material. 

The BR tag played a central role in the previous analysis M , since it was used for counting the number of 
prompt J/'ijj candidates. Subsequent improvements in ECAL stability and removal of material in the magnet 
region resulted in reduced combinatoric background and better mass resolution. Consequently the BR tag 
is now only used for calibration and cross-checks and not for signal counting. 

The invariant mass spectrum obtained for the final J/'ijj —>■ e+e“ sample is shown in Fig. As for 
the muon case, the number of prompt J/ijj candidates is determined by fitting the mass distribution with 
a function determined from MC studies. It includes 1) Gaussian functions for the J/'ijj and for the 'ijj{2S) 
signal, 2) a function to take into account the radiative tail of the charmonium decays and 3) an exponential¬ 
like function for the background. The radiative tail contains both the partially reconstructed J/ip 
decays as well as the tail due to incomplete electron energy measurement or BR recovery. The mean value 
of the J/ijj mass is 3109 ± Igtat MeV/c^, and the FWHM of the background subtracted J/'ijj peak is 140 
MeV/c^. The total number of J/'ijj e+e“ candidates is np(e+e“) = 103800 ± 1000. 

VI. DETACHED VERTEX ANALYSIS 

In the second stage of the analysis, J/'ijj mesons produced in the decay of a 6 hadron are searched for. 
The mean decay length of a 6 hadron in the HERA-B acceptance is about 8 mm. Accordingly, the J/'ijj 
mesons produced in b decays are typically well-detached from the primary interaction at the target wires. 
The b J/'ijj candidate selection is performed in two tightly correlated steps: first, a J/'ijj candidate (as 
described in the previous section) is assigned to one target wire and to one primary interaction vertex. Then, 
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FIG. 2: 



e^e invariant mass (GeV/c^) 

Invariant mass plots for both J/?/)-decay channels: a) muon sample, b) electron sample. 


if the association is unambiguous, cuts designed to distinguish J/ip mesons coming from b decays from those 
produced promptly at the target are applied. Almost all cuts are applied on the significance of the relevant 
quantity (ratio of a quantity to its uncertainty) rather than on its unnormalized value, to take advantage of 
the error information. The variables used to isolate detached J/ip mesons are described in the following. 

The decay length is defined as the distance between the primary interaction point and the J/ip vertex along 
the proton beam direction (which defines the HERA-B z axis). The positions of target wires, averaged over 
many consecutive events, define better the primary interaction point than the fitted positions of individual 
primary vertices; therefore the decay length is calculated with respect to the target wire. With a dilepton 
vertex resolution of about 450 /rm in z, the decays of b hadrons can be separated from primary vertices with 
high efficiency. 

For J/ip mesons produced at a detached vertex, each of the decay leptons is likely to be inconsistent with 
being produced at the target. Therefore a cut on their impact parameters to the wire, defined as the lepton 
distance from the wire at the z-coordinate of the wire, is applied. This requirement is very effective for 
reducing background, as the wire position is known in two dimensions with high precision. 

On the other hand, since the direction of the produced J/ip meson is nearly collinear with that of the b 
hadron, the impact parameter of the J/ip to the wire is typically small (mostly below 20/rm). An upper 
cut is applied to remove combinatorial background, as well as background from semileptonic bb and double 
semileptonic cc decays. 

For double wire runs and for multiple primary vertex events, the quantities described above can only be 
defined after association of the dilepton candidate to the correct target wire and to the correct primary 
interaction vertex. For this purpose, we define two x^-like quantities based on the lepton and the dilepton 
impact parameters to a wire and on their distances of closest approach to a primary vertex. The wire and 
the vertex for which the corresponding is minimum are associated to the dilepton candidate. Events in 
which the wire of the dilepton and the wire of the associated primary vertex do not match are discarded. In 
this way, the probability for wrong wire assignment, estimated from bb MC events, is less than 1%. 

Moreover, in the case of double wire runs, it is possible that a lepton candidate from an interaction on one 
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wire is combined with a lepton candidate originating from the other wire. Since the reconstructed vertex of 
these two leptons tends to lie between the two wires, these events are a dangerous source of background for 
the detached signal. To reduce this type of background, events in which one or both leptons are pointing to 
a reconstructed primary vertex on a different wire than the assigned one are discarded. 

Events having unambiguous wire-vertex association are analyzed further. We exclude events with the 
absolute value of the decay length |Az| < 2 mm (fiducial volume cut), in order to ensure that the J/V’ vertex 
is well outside of the target wire. This removes more than 99% of prompt J/4> background while keeping 
more than 75% of the signal. Afterwards, on the total sample a cut optimization procedure is applied, on 
the following quantities: the Az significance, the lepton and the dilepton impact parameters to the wire and 
the lepton and dilepton distances of closest approach to the primary vertex and E/p for electrons. 


A. Cut Optimization Procedure 

To avoid biasing our signal yield and resulting cross section ratio, we choose our selection criteria ’blindly’, 
i.e. without looking at events in the signal region until all selection criteria are finalized. The quantity to be 
maximized as a function of the cuts is the signal significance S/^/S BG in a fixed dilepton invariant mass 
window, where the signal (5) is taken from b J//} MC events, and the background {BG) is obtained by 
combining MC and real data. The main background sources are listed below. The first three contributions 
are estimated from MC, the fourth from data: 

• A possible prompt J/ip background surviving detachment cuts is one of our main concerns, since 
there is no way to distinguish it from J /ip mesons coming from b decays. Its yield is proportional to 
the prompt J/ip cross section. A large MC sample (« 10^ events) is used to estimate the surviving 
background contribution. 

• A large fraction of the background is coming from bb events where both b quarks produce a high-pr 
lepton. The yield of this background is proportional to the bb cross section. 

• Double semileptonic cc decays can contribute in a similar way. The intrinsic pr cut at the trigger level 
removes most of this background. Nevertheless this contribution cannot be neglected, since the cross 
section for cc events is larger than the bb cross section by three orders of magnitude. 

• Another contribution is due to combinatorial background, i.e. dileptons made of incorrectly recon¬ 
structed or fake tracks, or tracks coming from kaon or pion decays in flight. Its size can be estimated 
from the number of events in which the dilepton vertex is reconstructed in the unphysical region up¬ 
stream of the target. Also, like-sign lepton pairs (available only for the muon channel) or dileptons 
artificially formed by tracks from different events can give an estimate of this background. 

The background used for the cut optimization is an appropriately scaled combination of prompt J/ip 
MC, bb and cc double semileptonic MC events and either like-sign data events (in the muon channel) or 
upstream data events for the combinatorial background. The relative contribution of the different samples is 
estimated on the basis of the relative cross sections (for the MC samples) and from the upstream spectrum 
for the combinatorial background. Since the bb cross section is not known (being the goal of the current 
measurement), its value is left to vary in a range between a quarter and double of our previous measurement 
M- The optimization is repeated for several intermediate values. A stability check is also performed by 
varying the initial values of the cuts in order to verify the independence of the minimization results on the 
initial conditions. 

This approach can work effectively only if the data distributions are well reproduced by MC. This is 
verified for all the quantities used for the selections. As an example, in Fig. 13 the distributions of the lepton 
impact parameter to the wire are shown for a sample of background subtracted J/ip —>■ (a) and a 

sample of background subtracted double bremsstrahlung tagged J/ip ^ e+e“ (b). Agreement between data 
and MC is obtained for all variables used in the optimization procedure for both channels. 

Of the six parameters entering in the optimization procedure, the cuts on the lepton’s and dilepton’s 
closest distance of approach to the primary vertex are fixed by the optimization procedure to a value where 
all events surviving the other cuts are selected. Therefore these cuts have been removed from the optimization 
procedure and will not be mentioned further. 
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FIG. 3: Comparison between the distribution of the lepton impact parameter to the wire for real data (points 
with error bars) and for the prompt J/i/; MC (histogram), a) J/ip candidates selected in the region 

2.95 < m^+^- < 3.25 GeV/c^ after sideband subtraction for real data, b) double BR tagged » e'*"e“ candidates 

selected in the region 2.9 < m^+^- < 3.2 GeV/c^ after sideband subtraction for real data. The continuous solid line 
shows the b ^ J/ip MC sample (arbitrary normalization) while the arrows mark the cut position determined by the 
optimization procedure. 


B. J/ipX ^ n+^-X 

The cut values obtained with the optimization procedure for the muon channel are listed in the second 
column of Table Q] Once the cuts are applied, only 268 events with secondary vertices downstream of the 
primary interactions and having a dilepton invariant mass above 2.2 GeV j(? survive, while upstream of the 
target (unphysical region for b ^ J/4’ decays), 33 events are found in the same mass region. The invariant 
mass distribution of these detached candidates is shown in Fig. 01 where a clear peak corresponding to the 
J/ip mass is present only in the downstream sample. 


Cut 

fj. value 

e''"e value 

Optimized 

Absolute Az 

> 2 mm 

> 2 mm 

no 

Az significance 

> 9.0 a 

> 10.0 a 

yes 

Lepton impact to wire 

> 2.6 (7 

> 3.0 (7 

yes 

J/-!/! impact to wire 

< 9.0 (7 

< 12.0 CT 

yes 

Lepton identification 

> 0.05 

—2.0 a < Ejp — 0.99 < 3.5 a 

E/p only 


TABLE I: Cut values used to select the detached J/ij) candidates in the muon and in the electron channels. Quantities 
submitted to the blind optimization procedure are explicitly indicated in the last column. 

To count the detached J/'0 candidates, an unbinned maximum likelihood ht is performed using the recon¬ 
structed invariant mass values. In the fit, the function used for the prompt J/ij} analysis models the peak 
and an exponential function models the background. As free ht parameters, we use the J/ip and background 
yields and the slope of the exponential. When the mass position and width of the peak are left as free 
parameters, the obtained values are consistent with those of the prompt J/4> spectrum. Thus, the hnal 
results are given with the mass and width hxed to the values of the prompt signal. In Fig. the solid 
line shows the result of the likelihood ht. For the ht to the total data sample, we hnd = 46.2^® g and 
rihkg = 222 ±15 events for the signal and the background, respectively, and a background slope parameter 
of A = -1.9 ±0.1 (GeV/c^)-!. 

As a cross check, a ht to the upstream events is performed with the same free parameters. In the full 
sample, the ht hnds ng; = —0.1±1.4 J/'0 events, compatible with absence of signal, Ubkg = 33±6 background 
events and A = —1.9 ± 0.3 (GeV/c^)“^. 

Using the like-sign spectrum and the background contributions surviving the optimization procedure, the 
hnal composition of the background in the detached dilepton spectrum has been estimated. Within statistical 
uncertainty, all the background can be accounted for by the four background sources described above, with 
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FIG. 4: Dimuon mass distributions after vertex detachment cuts, from the full muon sample, a) Upstream events 
(combinatorial background), b) downstream events. The solid line shows the result of the likelihood fit. 


the dominant contributions being combinatorial (44%) and bb (43%) semileptonic decays. Since charm 
particles have a relatively short lifetime, their contribution to the background is suppressed (13%) by the 
detachment cuts. From MC, no prompt J/tp event (< 0.7 at 90% C.L.) is expected to survive the detached 
vertex selection. Finally, in the region upstream of the target the background is purely combinatorial. 

To determine Ra<t in our xp range according to Eq. the number of prompt J/ip mesons and the 
efficiency terms must be evaluated for each of the 14 different sub-samples (one sample for each wire in each 
wire configuration, labelled i later). The number of prompt J/ip —> candidates in each sub-sample is 

evaluated with the procedure described in Sect. EH To obtain the efficiencies, the MC events of prompt 
J/ip and b J/ip are subjected to the same analysis chain used for the data. The average final value of 
the efficiencies is (e_R,i • = 0.400 ± 0.006, where the uncertainty is given by the statistical fluctuations 

and sub-sample variability only. All the results obtained for different target materials are summarized in 
Table HU The resulting value for in our kinematic range, averaged over all sub-samples and target 

materials, is: 


= 0.0295 ± 0.0055stat (3) 

where the uncertainty is statistical only. The systematic uncertainties as well as the tests of the stability of 
the result are discussed in Sect. EnH 


C. J/ipX e+e-X 

As discussed in Sect. EH the J/tp ^ e+e decays are affected by a larger background and the selection 
of the detached vertex signal is therefore more critical. Since the electron particle identification (mainly 
based on the E/p ratio of the lepton tracks) is crucial for candidate selection, this quantity is included in 
the blind optimization procedure. The cut values obtained for the electron channel are listed in the third 
column of Tabled The Table shows that the optimized values for the cuts applied in the e+e” channel 
are close to those for the muon channel, albeit generally stronger. Only 229 events with dilepton vertices 
downstream of the primary interactions and with dilepton invariant masses above 2.0 GeV /(? survive these 
cuts, while 51 events with upstream dilepton vertices are found. The invariant mass distribution of these 
detached candidates is shown in Fig. [3 a peak is visible in the downstream sample at We+e- = 

The same type of unbinned maximum likelihood fit used in the muon channel is performed also on the 
detached electron candidates, and the fit result is shown in Fig. Eb. For the total data sample, we find 
rifjj = 36.9I7 I and ribkg = 192 ± 15 events for the signal and the background, respectively, and a background 
slope parameter of A = —1.7± 0.1 (GeV/c^)“^. 

As a further cross check, a fit to the upstream events is performed with the same free parameters. The fit 
finds = —4.2 ± 1.9 signal events (compatible with absence of signal events), Ubkg = 55 ± 8 background 
events and A = —1.6 ± 0.2 ( GeV /c^)~^. 
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FIG. 5: Dilepton mass distribution after vertex detachment cuts, from the full electron sample, a) Upstream events 
(combinatorial background), b) downstream events. The solid line shows the result of the likelihood ht. 


From the optimization procedure we can also extract the composition of the remaining background. As in 
the muon channel, the four sources considered can account for the obtained spectrum. The bb semileptonic 
decays are the dominant contribution (49%), followed by the combinatorial background (35%) and the open 
charm events (16%). Also, in this case no prompt J/?/’ event (< 1.1 at 90% C.L.) is expected to survive the 
applied cuts for b selection. 

As for the muon measurement, i?A(T is evaluated by determining the number of prompt J/ip mesons 
and the efficiency terms separately for each sub-sample (Eq. |2l). The average value of the efficiencies is 

j -eA^.) = 0.40 ±0.02, which is close to that obtained in the muon analysis. The resulting value for 
in the electron channel in our kinematic range is: 

RAa{e+e-) = = 0.0353 ± 0.0078,tat, (4) 

in good agreement with the muon result (see Table El). 


Channel 

fi+H 

e'’'e 

Target 

Carbon 

Titanium 

Tungsten 

Carbon 

Titanium 

Tungsten 

Atomic weight 

12.01 

47.87 

183.84 

12.01 

47.87 

183.84 

Prompt J/'ip (np) 

93700 ± 300 

8080 ± 100 

45560 ± 200 

67100 ± 700 

4800 ± 200 

32400 ± 600 

Detached J/'ip 

27.8 ±6.3 

3.0 ±2.1 

15.5 ±4.8 

17.81®:® 

0.9 ± 1.0 

18.41®:® 

• eft))) 

0.398 ± 0.004 

0.397 ±0.011 

0.404 ± 0.005 

0.366 ± 0.007 

0.424 ± 0.022 

0.402 ± 0.014 

a 


0.96 ±0.01 



0.96 ±0.01 


Br(66^ J/'ipX) 


2.32 ±0.20 



2.32 ±0.20 


(xlO-^) 

3.22 ± 0.73 

4.0 ±2.8 

3.6 ±1.1 

3.2 ±1.0 

1.8 ±2.2 

6.2 ± 2.0 

Aa{bb)/Aaj/^ (xl0“^) 


2.95 ±0.55 



3.53 ± 0.78 



TABLE II: Quantities entering into the cross section ratio measurements for the two channels and for the three 
target materials. All the numbers are given for the xf interval [-0.35,-1-0.15]. 


VII. FURTHER b FLAVOR CONFIRMATION ANALYSES 

To confirm the b flavor content of the selected detached vertex sample, several independent cross checks 
have been performed on it and on a second sample obtained by relaxing the selection cuts. The checks 
include an estimation of the mean lifetime of the detached candidates. Also a search is made for other decay 
products of the b hadrons in the enlarged selection, with the aim of a complete identification of exclusive 
decay channels or the identification of a potentially cleaner sample of three-prong detached vertices. These 
and other tests are described in the following. 
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proper time (ps) 



FIG. 6: Scatter plot of fj,^ ij.~ (a) and e"'"e“ (b) invariant mass versus the proper time for the detached candidates. 
The downstream detached candidates are marked with filled circles in the J/ij) regions (3.0 < m^+^- < 3.2 GeV/c^ 
and 2.9 < ni^+^- < 3.2 GeV/c^) and with open circles elsewhere. Upstream candidates are shown with open triangle 
marks. The hatches define the regions excluded from the mean lifetime measurement. 


A. Lifetime Fit 

The determination of the mean lifetime of the detached candidates is one of the clearest confirmations of 
the b flavor of the signal. A direct and precise measurement is however not possible, since in the inclusive 
6 —> J/'0 analysis the b hadron is not fully reconstructed and therefore its momentum and /Jy value are 
not known. However, MC studies show that a good estimator of the b hadron’s /Jy is the corresponding 
value of the J/'0 meson produced in the decay. Because of this strong correlation, the proper time ti of a 
b candidate can be approximated as ti = /S.Zi/((3^)j/^ ■ k, where Azi is the decay length, {l3j)j/.,p is the 
(3^ factor evaluated on the dilepton and fc is a numerical factor (« 1.03) which corrects for the decay and 
trigger bias. Within the limits of this approximation, a mean lifetime measurement is therefore possible even 
without exclusive identihcation of b hadrons. 

In Fig. 1^ the scatter plots of the dilepton invariant mass versus the proper time for the muon and the 
electron candidates are shown. As can be seen, the proper time distribution behaves differently depending 
on the invariant mass region. To quantify these differences, unbinned maximum likelihood hts in different 
regions have been performed, assuming the usual exponential decay distribution for the signal. In these hts, 
background events are treated like the signal: only one ht parameter, which represents the mean lifetime of 
the decay, is assumed in the likelihood function. The ht range is limited to the region where trigger and 
selection efficiencies are large: 0.9 - 8.5 ps. Since the efficiency is not constant in this range, its variation 
is included in the htting model, by weighting each event by the reciprocal of its detached vertex selection 
efficiency. The results of the hts are summarized in Table nni 

For detached candidates in the J/ip mass range ([3.0, 3.2] GeV/c^ in the muon channel and [2.9, 3.2] GeV/c^ 
in the electron channel), the mean lifetime is 1.61±0.27 ps and 1.21±0.18 ps for dimuons and dielectrons, 
respectively, in agreement with the value obtained from generated and reconstructed MG b decays: 1.56±0.01 
ps. Background regions in data (side bands in the invariant mass spectrum and upstream events) usually 
have lower mean lifetime values and larger uncertainties due to the poorer statistics. 

For reference we also measure the mean lifetime of the main background channels described by MC. 
As expected, the open cc events generally have a lower lifetime, while the semileptonic bb events behave 
differently depending on the invariant mass range of the dilepton. The higher mass range is dominated by 
combinations of leptons from two diherent b decay branches, while, in the low mass range, dileptons from 
the decay chain of a single b quark dominate. The mean lifetime measured in the latter sample is clearly 
close to that of the b. 

To improve the statistical accuracy of the lifetime measurement, a joint unbinned likelihood fit is performed 
on the muon and electron candidates (104 events in total). This sample is highly enriched in b content. 
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Sample 

Mass Range (GeV/c^) 

(ps) 

Tee (ps) 

data downstream 

J/ip region 

1.61 ±0.27 

1.21 ±0.18 

b J/ip MC 

J/ip region 

1.55 ±0.01 

1.57 ±0.01 

data downstream 

2.0-2.5 

1.30 ±0.16 

0.86 ± 0.09 

data downstream 

3.6-12.0 

0.55 ±0.13 

0.54 ±0.21 

data upstream 

2.0-12.0 

1.08 ±0.21 

1.16 ±0.19 

bb semileptonic MC 

2.0-2.5 

1.47 ±0.12 

1.45 ±0.04 

bb semileptonic MC 

3.6-12.0 

0.71 ±0.09 

0.67 ±0.03 

cc double semileptonic MC 

2.0-12.0 

0.47 ±0.06 

0.30 ± 0.02 


TABLE III: Results of the unbinned maximum likelihood fits on the detached j-i and e''"e candidates in real 
data and on various MC samples. The J/ijj mass range is defined as for Fig. 



FIG. 7: a) Likelihoods for the and e''"e detached candidate fits and joint likelihood as a function of the mean 
lifetime; b) distribution of proper times for the detached J/'i/) events corrected for selection efficiency. 


having about 80 6 ^ J/ij) candidates plus about 10 events coming from semileptonic b decays. A value of 

= 1.41 ± 0.16 ps is obtained (see Fig.[7^). Since the total number of candidates is sufficiently large, a 

standard fit is also performed on the efficiency-corrected distribution of proper times for the same events 

2 

(see Fig. m- In this case, the minimum (y^/ndf = 5.1/8) is obtained for = 1.39 ± 0.19 ps. Both 
results are in good agreement with the expected value for b hadrons. 


B. J/ip + analysis 

As described in Sect. EH rather strong (and therefore relatively inefficient) detachment cuts are necessary 
to do an inclusive selection of 6 ^ J/ip decays. A partially different sample of b events can be obtained by 
relaxing the detachment cuts and imposing additional requirements on other quantities which characterize 
b hadron decays, for example requiring the presence of a third track forming a good vertex with the 
J/'0 candidate. Since the J/ip carries most of the b hadron momentum, the additional track usually has 
low momentum and points far away from the primary interaction. This selection gives access to b events 
with smaller detachment from the target wire, but nevertheless provides high purity because of the stringent 
requirements on the three-prong decay vertex (which has an intrinsically lower background). 

For events with dimuon candidates selected as described in Sect. m a search is made for additional decay 
products among good quality tracks nearby the J/4> vertex. The presence of a third track, whose closest 
distance of approach to the primary vertex is larger than 300 /rm is required. Detachment cuts are then 
applied, and the requirements summarized in Tableware relaxed. In particular the minimum Az significance 
is reduced to her, and the minimum lepton impact parameter to the target wire is lowered to 2a. 

The dilepton invariant mass and mass vs proper time distribution of events selected in the muon sample 
are shown in Fig.|3 A J/ip peak is clearly visible, with an improved signal-to-background ratio with respect 
to the detached J/ip analysis of Sect. I VII Fig.E From the unbinned maximum likelihood fit, the number of 
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invariant mass (GeV/c^) jiV invariant mass (GeV/c") 



FIG. 8: invariant mass for the J/tp + analysis for upstream (a) and downstream (b) events. The line in 

(b) shows the result of the unbinned max likelihood fit. c) invariant mass as a function of the proper time for 

the detached candidates. 


detached J/ip events is nj^i = 22 ± 5. Only 50% of these events are in common with the previous detached 
vertex analysis. 

Using these selected events, another value of the cross section ratio in our acceptance, TJact = 
Aa{bb)/Aaj/^ = 0.043 ± 0.010, was obtained. The measurement is compatible with both the muon and 
the electron results presented in Sect. I VII 1 within 0.6 ct from the electron value, obtained with a statistically 
independent sample). Since it is partly correlated with the muon result of Sect. IVl'El it is not combined into 
the final cross section results presented in Sect. Emi 


C. other Studies 

Further confirmation of the b flavor content of the detached J/ip sample has been obtained from other 
analyses. We studied the kaon population and the distribution of high pT tracks. For muon events and for 
two-BR tag electron events falling in an invariant mass window around the J/ip ■, we observed that the real 
data after prompt selection behave as the prompt J/ip MC candidates, while after the detachment selection 
larger values of both, the pT distribution and the kaon population, are found, compatible with the 5 —> J/'0 
MC sample. 

A search for the exclusive decays —> J/ipK~^ and J/ipK'^TT~ has been made. A few fully 

reconstructed B meson candidates are found, both for the muon and the electron channels, but not sufficient 
for a full exclusive analysis. 

The kinematic characteristics of the detached J/ip candidates have been checked against the production 
model employed to evaluate the efficiencies and were found to be compatible in our xp and pr acceptance. 


VIII. COMBINED CROSS SECTION MEASUREMENT 

The two measurements of Aa{bb)/Aaj/^, presented in Eqs. © and 0 , are statistically independent 
and compatible. To have a more precise measurement of Rao-i a joint unbinned maximum likelihood fit is 
performed on the detached pAp.~ and e+e“ candidates, using Rao- and all background terms (ribkg and A 
for the two channels) as free parameters. In Fig.|^ the likelihoods of the invariant mass fits are shown for 
pApL~ and e“'"e“ separately and for the combined fit. The fit provides the final value 

RAa = = 0.0314 ± 0.0049,tat 
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FIG. 9: Likelihoods of the invariant mass fits as a function of i?Acr, shown separately for data (dotted line), 

e^e“ data (dashed line) and for the combination of the two (continuous line). The dot-dashed line shows the likelihood 
function corresponding to the 2000 data analysis which is used in the final average. 


where the quoted uncertainty (15%) is dominated by the statistical fluctuations on the detached J/'ijj counting 
and contains the statistical contributions from the prompt J/ijj counting (f« 1%), the efficiencies (~ 5%) and 
the a exponent (« 1%). 

The main sources of systematic uncertainties^, which are not related to the 66 statistics, are due to the 
Br(66—> JlijjX) (8.6%), to the trigger and reconstruction efficiency ratio er (5%), to the 6 production and 
decay model (5%), to the prompt J/^j} production model (3.1%) and to the prompt J/iji counting (1.5%). 
Other contributions are below the 1% level. To determine the sensitivity of the results to the cut values, we 
vary them within reasonable ranges, always requiring a negligible prompt J/ip background. A conservative 
estimate of 5% systematic uncertainty has been obtained. The uncertainties in the background shape and 
yield give a negligible contribution to the systematic uncertainty (below 1%) in the muon channel but a 
sizeable contribution in the electron channel (7%). The overall systematic uncertainty for the two channel 
average measurement is 14%. 

The weighted averages of the muon and electron results for the different target materials and for the total 
are presented in Table ESI As can be seen from the probability, the two channels are always compatible. 




The three measurements are also compatible with the assumed linear A dependence of the 66 cross 


section. 

The present result is l.TcTstat below the previous HERA-B measurement [T3 (Ra. = 0.107 ±0.047stat± 
0.022sys) which was obtained with a similar technique but with a much smaller data sample (dot-dashed line 
in Fig. 0. Since the two data samples are independent, we combine the two measurements assuming that 
the systematic errors are uncorrelated. We obtain: 


77 a CT 


Aa{bb) 

Aajj^ 


0.032 ± 0.005stat ± 0.004sys. 


(5) 


Material (A) 

Vj/v. 

'^bb 

a ratio value 

Prob(x^) 

Carbon (12) 

161 k 

45.6 ±8.4 

0.0324 ± 0.0059 

94% 

Titanium (48) 

13 k 

3.9 ±2.3 

ActA/Act^ 0.027 ± 0.018 

54% 

Tungsten (184) 

78 k 

33.9 ± 7.6 

0.044 ± 0.010 

29% 

Total 

251k 

83± 12 

Aa{bb)/Aaj/^ 0.0314 ± 0.0049 

54% 


^ The quoted systematic uncertainties correspond to a Icr interval and are usually evaluated as the maximal interval width of 
the cross section variation divided by \/T2. 
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TABLE IV: Combined j-i and e"'"e results: sum of the number of prompt {Nj/^) and detached J/’iji found 
for the two channels together, Aa^/Aap (rows 1-3 ), Aa(bb)/A gj/^ (row 4), and probability of the average 
value. 


To compare our result with other measurements and theoretical predictions, we extrapolate the i?A(T ratio 
to the full kinematic range and then make use of the value of presented in Sect.O obtaining the cr(65) 
value: 


aibb) = 14.9 ± 2.2stat ± 2.4sys nb/nucleon. (6) 

This result can be compared with the available measurements, which were obtained in pAu Q and pSi Q 
collisions at 800 GeV/c proton momentum. Due to the near-threshold conditions in which all experiments 
(including HERA-B) operated, it is necessary to rescale the Fermilab results for the different \/s value of the 
collision. The latest QCD calculations medict an increase of the bb cross section between 800 (Fermilab) and 
920 GeV/c (HERA-B) by (42±9)% @,l3- All the available measurements are presented in Table IVl where a 
scaling of the Fermilab measurements to HERA-B energies is performed. As can be seen, the present result 
(Eq. inil is consistent with both the E789 value (within l.Gcr) and with the E771 value (within 1.8 (t). 


Expe¬ 

riment 

Year 

Target 

V 

(GeV/c) 

Events 

a{bb) 

(nb/nucleon) 

a{bb) 

at 920 GeV/c 

Ref. 

E789 

E771 

HERA-B 

HERA-B 

1995 

1999 

2002 

2005 

Au 

Si 

C/Ti 

C/Ti/W 

800 

800 

920 

920 

19 

15 

10 

83 

5.7 ±1.5 ±1.3 
43ti/ ± 7 

00+14 +6 (*) 

'^^-12 -7 

14.4 ±2.2 ±2.3 

8.1 ±2.2 ±1.9 

61/^1 ± 11 
00+14 +6 (.) 
' 44:_12 -7 

14.4 ±2.2 ±2.3 

m 

t.w. 

HERA-B 

2002/5 

C/Ti/W 

920 

93 

14.9 ± 2.2 ± 2.4 

14.9 ±2.2 ±2.4 

t.w. 


TABLE V: Present experimental situation on the bb cross section in pN interactions with the comparison of all the 
measurements rescaled at 920 GeV/c proton momentum. 

Value as quoted in the paper, which was obtained using a aj/^ value different from the current paper. 


The measurements presented in Table |V| can also be compared to the latest QCD predictions HQ , which 
have been performed by two theory groups. They both employ Next-to-Leading-Order (NLO) calculations 
with resummation techniques to take into account the large corrections due to emission of soft gluons. The 
calculations give a{bb) = 28 ± 15 nb/nucleon u and a{bb) = 251^3 nb/nucleon ^3|' the results are shown 
in Fig. cni as a function of ^/s. As can be seen, the present measurement is compatible with the theoretical 
calculations. 


IX. CONCLUSIONS 

A search for b l^l~X decays has been performed on a sample of 164 million dilepton trigger 

events, acquired during the HERA-B 2002-2003 physics run. The data analysis, based on identification of 
detached vertices, resulted in 46.2 ±8.2 b ^ J/'ipX candidates and 36.9 ±8.1 b ^ J/%I)X 

e+e“X candidates. 

From these events, we measure the ratio Aa{bb)/Aaj/^ of the bb production cross section to the prompt 
J/ip cross section in the HERA-B acceptance range (—0.35 < xp < 0.15). Within statistical uncertain¬ 
ties, the results for the muon and for the electron channel are compatible and yield a combined value of 
Aa{bb)/Aaj/^ = 0.0314 ± 0.0049stat ± 0.0044sys, where the systematic uncertainty (14%) is dominated by 
the 8.6% contribution due to the branching fraction Bi{bb^ J/ipX). The present measurement is then 
combined with the HERA-B 2000 measurement ^3: for a final result of: 

i?Aa = = 0.032 ± 0.005stat ± 0.004sys. (7) 

Aojj^ 

Further investigations of the detached vertex sample confirm the b nature of the candidates. The measured 
mean lifetime r = 1.41 ± 0.16 ps agrees well with previous measurements of the average b lifetime m. A 
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Proton energy (GeV) 


FIG. 10: Comparison of the available a{bb) measurements with the theoretical predictions of R. Bonciani et al. 
updated with the NNLO parton distribution function in 1^ (solid line: central value, dashed lines: upper and lower 
bounds) and N. Kidonakis et al. 0 (dot-dashed line: central value, dotted lines: upper and lower bounds). The 
upper and lower bounds have been defined by changing the renormalization and factorization scales (by a factor 0.5 
and 2.0) and the b quark mass (4.5 — 5.0 GeV/c^). 


search for extra tracks coming from the detached J/ip vertex in the muon channel yields 22 ± 5 partially 
reconstructed b^J/%j^ + K^-\-X events. From this latter sample, 50% of which consists of detached J/tp 
mesons not found in the inclusive analysis, a i?Acr cross section ratio in good agreement with the main 
measurement is obtained. In addition to these checks, we have verified that the detached candidate events 
have a kaon and pt distribution that is in good agreement with that expected for b events. 

To compare our result to previous measurements and to theoretical calculations, we compute the a{bb) 
cross section by extrapolating the ratio i?Acr to the full phase space and using the value of = (502 ± 
44) nb/nucleon ^3. The resulting <j{bb) cross section is 14.9±2.2stat±2.4sys nb/nucleon. Comparing to the 
other available experimental results, the present value is within 1.6cr of the E789 value Q (after rescaling to 
the same -^/s) and l.Str below the rescaled E771 measurement [ij. 

The HERA-B value can be compared with the latest QCD predictions, calculated beyond the NLO order 
and considering the effects of soft-gluon re-summation, performed at 920 GeV/c proton momentum. Our 
result is within Icr of the current calculations 0,113 
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